INTEGRATION OF A FINITE ELEMENT MODEL WITH THE DAP BONE REMODELING MODEL TO 
CHARACTERIZE BONE RESPONSE TO SKELETAL LOADING 

C.R. Werner 1 , L. Mulugeta 2 , J.G. Myers 3 , and J.A. Pennline 3 

! ZIN Technologies, 6745 Engle Road, Airport Executive Park, Cleveland, OH 44130 
2 Universities Space Research Assoc., Div. of Space Life Sciences, 3600 Bay Area Blvd., Houston, TX 77058 

3 NASA Glenn Research Center, Cleveland, Ohio 

INTRODUCTION: NASA’s Digital Astronaut Project (DAP) has developed a bone remodeling model that has been 
validated for predicting volumetric bone mineral density (vBMD) changes of trabecular and cortical bone in the 
absence of mechanical loading [1]. The model was recently updated to include skeletal loading from exercise and free- 
living activities to maintain healthy bone using a new daily load stimulus (DLS). This new formula was developed 
based on an extensive review of existing DLS formulas [2], as discussed in the abstract by Pennline et al. The DLS 
formula incorporated into the bone remodeling model utilizes strains and stress calculated from finite element model 
(FEM) of the bone region of interest. The proximal femur was selected for the initial application of the DLS formula, 
with a specific focus on the femoral neck. 

METHODS: The FEM was generated from CAD geometry of a femur using de-identified CT data [3, 4]. The femur 
was meshed using linear tetrahedral elements Figure (1) with higher mesh densities in the femoral neck region, which 
is the primary region of interest for the initial application of the DLS formula in concert with the DAP bone remodeling 
model. Nodal loads were applied to the femoral head and the greater trochanter and the base of the femur was held 
fixed. An L2 norm study was conducted to reduce the length of the femoral shaft without significantly impacting the 
stresses in the femoral neck. The material properties of the FEM of the proximal femur were separated between 
cortical and trabecular regions to work with the bone remodeling model. Determining the elements with cortical 
material properties in the FEM was based off of publicly available CT hip scans [4] that were segmented, cleaned, 
and overlaid onto the FEM. 

The FEM was solved using MSC Nastran 
101 linear static solution, and the output 
stresses were provided to the bone 
remodeling model to consolidate and 
determine a DLS outcome. The system of 
passing information from the FEM to the 
bone remodeling model was automated 
using Matlab. With the bone densities 
provided by the user, the FEM is initially 
read into Matlab and the initial material 
properties are written into the FEM file and 
solved using Nastran. Once complete, the 
output stresses are automatically read into 
Matlab and the DLS is calculated for a given 
exercise applied to the bone. As the bone 
density and bone volume fraction changes 
depending on overload or disuse, the FEM Figure 1: FEM of proximal femur (left) and block diagram of Matlab - 
is updated with new material properties Nastran interaction for the bone remodeling model (right) 

(Figure 1), which are functions of ash 

density and bone volume fraction, and the cycle is repeated until the duration of the simulation is complete. 
CONCLUSIONS AND FUTURE WORK: A FEM was successfully integrated with the DAP bone remodeling 
model DLS formula to simulate the influence of mechanical stimulus on the bone remodeling process in a fully 
automated manner. However, there are several improvements that will be made to increase the accuracy of the FEM 
such as using CT data to directly map the vBMD on an elemental level and applying anisotropic material moduli 
instead of the current assumption of isotropic material properties. Also, the remodeling algorithm and DLS formula 
may be applied to each element to evaluate vBMD changes per-element resolution rather than an aggregate of the 
trabecular and cortical regions. 

REFERENCES: [1] Pennline J and Mulugeta L (2014), 44th ICES , ICES2014-083; [2] Pennline J and Mulugeta L 
(2014), NASA Tech Memo , NASA/TM -20 14-2 18306; [3]; Werner C and Gorla RSR (2013), Int. J. Appl. Mech. Eng., 
18:911-921. [4] Harris MD et al. (2012) JOrthopRes, 30:1133-9. 
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I NASA’s Digital Astronaut Project Vision Statement 

I The Digital Astronaut Project (DAP) implements well-vetted computational models to predict and a 
I performance risks, and enhance countermeasure development. 

HRP Risks Knowledge Gaps Addressed by This Work 

1 Osteo4: We don’t know the contribution of each risk factor on bone loss and recovery of bone strength and which factors are the 
best targets for countermeasure application 

1 Osteo7: We need to identify options for mitigation of early onset osteoporosis before, during, and after spaceflight. 

1 Gap Fracture 3: We need a validated method to estimate the Risk of Fracture by evaluating the ratio of applied loads to bone 
fracture loads for expected mechanically-loaded activities during and after a mission 


OBJECTIVES 


I Model and Simulation Description 

1 The DAP has developed a bone remodeling model that has been validated for predicting 
volumetric bone mineral density (vBMD) changes of trabecular and cortical bone in the 
absence of mechanical loading [1] 

1 The model was recently updated to include skeletal loading from exercise and free-living 
activities to maintain bone using a new daily load stimulus (DLS) 

1 The DLS formula incorporated into the bone remodeling model utilizes strains and stress 
calculated from a finite element model (FEM) of the bone region of interest 

1 The proximal femur was selected for the initial application of the DLS formula, with a 
specific focus on the femoral neck 

| Objective 

' The goal of using the FEM in concert with the bone remodeling model is to broaden the 
capability of the bone remodeling model to included both loading and non-loading 
scenarios on bone. 

• The initial objective of combining the FEM with the bone remodeling model is to create a 
robust and fully automated simulation that helps serve as a research tool for bone. 


I FEM Description 

' The FEM is created based on the CAD of an anonymous subject’s CT hip scan [2] 

' Given the bone remodeling model is focused on the stress/strain around the femoral 
neck, the FEM of the CAD model is shortened in the long bone portion of the femur 
1 An L2 Norm study is conducted to ensure the stresses around the femoral neck are not 
affected by reducing length 

' The mesh is sliced into 1 inch sections and femoral neck stresses are used for 
comparison between iterations 


I FEM of Proximal Femur I 




METHODS 
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1 The volume is meshed using Siemens Femap meshing software (Munich, Germany) 

1 Element sizes were reduced around the femoral neck to provide better resolution at the 
location where the bone remodeling model is primarily focused 
1 The current approach of the FEM: 

■ Linear isotropic modulus of elasticity 

■ Linear tetrahedral elements 

■ NASTRAN 101 linear static solution 

■ Femap as the FEM pre and Post Processor 




Cortical Thickness Determination 



1 The selected cortical region is then cleaned up using 
Meshlab software, and imported into Femap in the form of 
null planar elements 

1 Since the FEM and cortical surfaces come from two 
different patients, it is assumed they will not perfectly align 
so the cortical thickness is scaled, oriented, and aligned to 
the FEM 

1 The FEM is then manually separated using freehand 
selection into cortical and trabecular regions using the 
cortical thickness model as a guide 


Loads and Boundary Conditions 

• The loading conditions are currently 
focused on mechanical stimulus due to 
walking and running 

• Two forces are applied: joint contact force 
at the femoral head and muscle forces at 
the greater trochanter 

• Used Orthoload data for joint contact 
load at the femoral head [3] 

• Muscle forces were estimated based 
on the method of Hazelwood and 
Castillo [4] 

• To get a better representation for the 
cumulative loading for a single load cycle, 
the mean force on the hip is used 

• FEM is constrained at the base. 


Daily Load Stimulus (DLS) 

k 

DLS £ S(<jj, N p fj , Pj , rj (Dj ), (Oj) . 
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• The DLS is used to convert the stress, as 
well as other parameters, to a singular useful 
value for the bone remodeling model 

• To develop this new formula, we conducted 
an extensive study to evaluate existing DLS 
formulas [5] 

Bone Remodeling Integration 

• The FEM is integrated by: 

• Receiving the modulus and load values 
from the remodeling algorithm 

• Determining the resultant stresses / 
strains 

• Passing them back to the remodeling 
simulation to calculate mechanical 
stimulus 

Nastran can be run multiple times 
throughout the simulation to update the 
load path due to changes in cortical and 
trabecular stiffness 
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F x_gt - C fgt * F x 
F y.gt = c fgt * F y 


F z_gt ~ Cfgt * F z 
F r_gt = C fgt * F 
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m Hazelwood and Castillo [4] 


k = The number of loading conditions/activities, or the number of exercises per day 
Nj = The number of loading cycles per loading condition or repetitions per exercise 
= Maximum principle stress per loading condition 
fj = The load frequency of a single repetition not including pause durations between reps. 
Pi = Repetition period, including any pause periods in between repetitions 
t j = The ratio between the total duration of the exercise or activity (Dj), over a full day 
Dj = Total duration of the exercise or activity bout 
o)j = Ratio of the current exercise speed over the nominal exercise speed 
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Reed FEM Input file and ImitlaHre 
FEM with Material Properties 



RrmodHInJ 
Model With 


A 


| 


1 1 
1 

Cortical bone 


"I 

1 

1 



MODEL VERIFICATION AND VALIDATION 


Analysis on change in bone density due to normal ** 
daily walking for different body masses and steps 
per day 

Results show no change in density with small * 
decreases for minimal steps per day cases 


Validation 

• Assumed that astronauts are likely to engage in 
higher daily activity that would be the equivalent 
to 12,500 to 18,000 steps 

• Mean body weight of 725 N and walking speed 
of 5km/h was assumed 

• The simulation results match the data reported 
by Lang et al [6] within the experimental 
standard deviation and standard error 

■ Post flight QCT data from 1 6 crewmembers 
collected at R0 and R =12 months reported 
in Lang et al [6] 
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• The cortical thickness determination of the bone required 
publically available CT scan data and segmented using 
3D slicer software (above) 

• The cortical region is selected using a combination of 
thresholding and manual selection to segment the cortical 
region in the bone 

• Resulting selection is exported as a raw .stl file 
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